A simple mechanism for the reversals of Earth's magnetic field 
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We show that a model, recently used to describe all the dynamical regimes of the magnetic field 
generated by the dynamo effect in the VKS experiment [f[ , also provides a simple explanation of the 
reversals of Earth's magnetic field, despite strong differences between both systems. 
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The Earth's magnetic field can be roughly described as 
a strong axial dipole when averaged on a few thousands 
years. As shown by paleomagnetic records, it has fre- 
quently reversed its polarity on geological time scales. 
Field reversals have also been reported in several nu- 
merical simulations of the geodynamo 0] and more re- 
cently, in a laboratory experiment involving a von Kar- 
man swirling flow of liquid sodium (VKS [H). It is worth 
pointing out that numerical simulations are performed in 
a parameter range orders of magnitude away from real- 
istic values, and that both the parameter range and the 
symmetries of the flow in the VKS experiment strongly 
differ from the ones of the Earth's core. We thus ex- 
pect that if a general mechanism for field reversals ex- 
ists, it should not depend on details of the velocity field. 
This is expected in the vicinity of the dynamo threshold 
where nonlinear equations govern the amplitudes of the 
unstable magnetic modes. We assume that two modes 
have comparable thresholds. This has been observed for 
dipolar and quadrupolar dynamo modes Q and has been 
used to model the dynamics of the magnetic fields of the 
Earth or the Sun [5|. However, in contrast to these pre- 
vious models, we consider two axisymmetric stationary 
modes and expand the magnetic field B(r, t) as 

B(r, t) = o(t)B x (r) + 6(t)B 2 (r) + . . . . (1) 

We define A(t) = a + ib and write the evolution equa- 
tion for A using the symmetry constraint provided by 
the invariance B — > — B of the equations of magnetohy- 
dro dynamics. This imposes A — > —A, thus the amplitude 
equation for A is to leading nonlinear order 

A = fiA + vA + 7i A 3 + -f 2 A 2 A + -f 3 AA 2 + 74 A 3 (2) 

where /1, v and 7^ are complex coefficients. Equations 
of the form ([2]) arise in different contexts, for instance 
for strong resonances, and their bifurcation diagrams are 
well documented [6[. Defining A = RexpiO, a further 
simplification can be made when the amplitude R has a 
shorter time scale than the phase 9 and can be adiabati- 
cally eliminated. In that case, 9 obeys an equation of the 
form 

9 = a a + (a„ cos 2n0 + f3 n sin 2n0) . (3) 

n>l 




FIG. I: Phase space of a system invariant under B — » B 
and displaying a saddle-node bifurcation: (a) below the on- 
set of the saddle-node bifurcation, (square blue): stable fixed 
points; (red circle): unstable fixed points, (b) Above the 
threshold of the bifurcation, the fixed points (empty symbols) 
have collided and disappeared, the solution describes a limit 
cycle. Note that in (a), below the onset of the saddle-node 
bifurcation, fluctuations can drive the system from Bs to Bu 
(phase Bs — > Bu) and initiate a reversal (phase Bu — * —Bs) 
or an excursion (phase Bu —> Bs) . 



The absence of odd Fourier terms results from the invari- 
ance B — > — B that implies 9 — > 9 + it. Stationary solu- 
tions of ^ disappear by saddle-node bifurcations when 
parameters are varied. When no stationary solution ex- 
ist any more, a limit cycle is generated which connects 
the former stable point 9 S to 9 S + ir, i.e., B s to — B s (see 
Fig. 1). This elementary mechanism for reversals is not 
restricted to the validity of but results from the two 
dimensional phase space of d2J) Thus, the qualitative 
features of the dynamics can be captured using the sim- 
plest possible model keeping the leading order Fourier 
coefficients ao and Pi (a\ can be eliminated by changing 
the origin 9 -> 9 + 9 ). 

So far, we did not consider possible effects of fluctua- 
tions. The flow in the Earth's core, as well as in the VKS 
experiment, is far from being laminar. We can therefore 
assume that turbulent fluctuations act as noisy terms in 
the low dimensional system that describes the coupling 
between the two magnetic modes. In Fig. la, the system 
is below the threshold of the saddle-node bifurcation and 
in the absence of fluctuation exhibits two stable (mixed) 
solutions. If the solution is initially located close to one 
of the stable fixed points, say Bs, fluctuations can push 
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the system away from Bs. If it goes beyond the unsta- 
ble fixed point Bu, it is attracted by the opposite fixed 
point —Bs, and thus achieves a polarity reversal. A re- 
versal is made of two successive phases. The first phase 
Bs — > Bu in Fig. la is the approach toward an unsta- 
ble fixed point. The deterministic dynamics acts against 
the evolution and this phase is slow. The second phase 
Bu — > —Bs, is fast since the deterministic dynamics fa- 
vors the motion. 

At the end of the first phase, the system may return 
toward the initial stable fixed point (phase Bu — > Bs), 
which corresponds to an excursion. We emphasize that, 
close enough to the saddle-node bifurcation, reversals re- 
quire vanishingly small fluctuations. To take them into 
account, we modify the equation for 9 into 



a a + ai sin(20) + A((t) , 



(4) 



from which we derive the evolution of the dipole by 
d = i?cos(f? + 9q). C is a Gaussian white noise and A 
is its amplitude. We have computed a time series of the 
dipole amplitude for a system below the threshold of the 
bifurcation {a x = -185 Myr" 1 , a /ai = -0.9, 6 = 0.3) 
and with noise amplitude A/yjoi | = 0.2. Note that a,\ 
is arbitrary at this stage. Its value results from a fit of 
paleomagnetic data (see below). The dipole amplitude 
is displayed in Fig. 2 together with a time series of the 
magnetic field measured in the VKS experiment and the 
composite record of the geomagnetic dipole for the past 
2 Myr. The three curves display very similar behaviors 
with abrupt reversals and large fluctuations. We have 
checked that similar dynamics are obtained when equa- 
tion ([2]) with noisy coefficients is numerically integrated. 

One of the most noticeable features common to these 
three curves is the existence of a significant overshoot 
that immediately follows the reversals. In Fig. 3, the 
enlarged views of the period surrounding reversals and 
excursions also show that this is not the case for the 
excursions. In fact, the relative position of the stable and 
unstable fixed points (Fig. 1) controls the evolution of 
the field. During the first phase, reversals and excursions 
are similar, but they differ during the second phase. The 
synopsis shows that the reversals reach the opposite fixed 
point from a larger value and thus display an overshoot 
while excursions do not. 

Below the onset of bifurcation, reversals occur very 
seldomly, which indicates that their occurence requires 
rarely cooperative fluctuations. The evolution from the 
stable to the unstable fixed point (phase Bs — > Bu in 
Fig. 1) can be described as the noise driven escape of the 
system from a metastable potential well. The durations 
of polarity intervals are equivalent to the exit time and 
are exponentially distributed 0] according to 



P[T] ocexp(-T/(T)). 



(5) 



The averaged duration (T) depends on the intensity of 
the fluctuations and on the distance to the saddle-node 
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FIG. 2: (Top) Time series of the dipole amplitude for a system 
below the threshold of the bifurcation (see text for the values 
of the parameters) (Middle) Time series of the magnetic field 
measured in the VKS experiment for the two impellers ro- 
tating with different frequencies Fl=22 Hz, F2=16 Hz (data 
from 0]). (Bottom) Composite paleointensity curve for the 
past 2 millions years, present corresponds to t=0 (data from 



bifurcation. For the model studied, we obtain 
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(6) 



which corresponds to (T) ~ 170 kyr for the parame- 
ters used in Fig. 2. We observe that deterministic pa- 
rameters are of the order of the Ohmic dissipation time 
(7r/|ai| ~ 17000 years) whereas much larger time scales 
are measured for (T) because of the low noise inten- 
sity. This explains that the mean duration of phases with 
given polarity is much larger than the one of a reversal. 
The above predictions assume that the noise intensity 
and the deterministic dynamics do not vary in time. It is 
likely that the Rayleigh number in the core and the effi- 
ciency of coupling processes between the magnetic modes 
have evolved throughout the Earth's history. The expo- 
nential dependence of the mean polarity duration (T) on 
noise intensity implies that a moderate change in convec- 
tion can result in a very large change of (T) . This might 
account for changes in the rate of geomagnetic reversals 
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FIG. 3: Comparison of reversals and excursions in the solutions of the equation for 9 (left), the VKS experiment (middle) (data 
from 0]), and in paleomagnetic data (right) (data from §]). Black curves represent the averaged curve, each realization being 
represented in grey. 



and very long periods without reversals (so called super- 
chrons). The reversal rate reported in Q is displayed in 
Fig. [?] together with a fit using Eq. [5] assuming a linear 
variation in time of the coefficients governing the distance 
between Bs and Bu. A simple variation of one parame- 
ter captures the temporal evolution of the reversal rate. 
Thus, although it can be claimed that there are several 
fitting parameters, this quantitative agreement strength- 
ens the validity of our model of reversals. 
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FIG. 4: Variation of the reversal rate (T) 1 close to a super- 
chron together with the fit with equation ([6]). Data (•) have 
been extracted from Q. |ai| = 185 Myr" 1 , A/,/ja7i = 0.2, 
ao/ai = -0.9(1 - r/650) for r < 80, a /ai = -0.9(0.55 + 
t/360) for t > 120, where r is time in millions of years. 

We now consider the different possibilities for the mode 
Ba(r) coupled to the Earth's dipolar field Bi(r). As- 
suming that the equator is a plane of mirror symmetry, 



the different modes can be classified as follows: dipo- 
lar modes are the ones unchanged by mirror symme- 
try, D — > D, whereas quadrupolar modes change sign, 
Q — > Q. From an analysis of paleomagnetic data, Mc- 
Fadden et al. have proposed that reversals involve an 
interaction between dipolar and quadrupolar modes [1 01 ] . 
In that case, Bi(r) and B2(r) change differently by mir- 
ror symmetry. If the flow is mirror symmetric, this im- 
plies that Eq. should be invariant under A — > A 
which amounts to 9 — > —6. Consequently, a n = in ([3]) 
and no limit cycle can be generated. We thus obtain an 
interesting prediction in that case: if reversals involve a 
coupling of the Earth's dipole with a quadrupolar mode, 
then this requires that the flow in the core has broken 
mirror symmetry. This mechanism explains several ob- 
servations made in numerical simulations: reversals of 
the axial dipole, simultaneous with the increase of the ax- 
ial quadrupole, have been found when the North-South 
symmetry of the convective flow is broken It has 

been shown that if the flow or the magnetic field is forced 
to remain equatorially symmetric, then reversals do not 
occur [l3| . The possible effects of heterogeneous heat flux 
at the core-mantle boundary (CMB) on the dynamics of 
the Earth's magnetic field have also been investigated 
numerically [T3|. Compared to the homogeneous heat 
flux, patterns of antisymmetric heterogeneous heat flux 
were shown to yield more frequent reversals. Within our 
description, this appears as a direct consequence of the 
breaking of hydrodynamic equatorial symmetry driven 
by the thermal boundary conditions. From the point 
of view of the observations, little is known on the ac- 
tual flow inside the Earth's core. It has recently been 
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noted that the ends of superchrons are followed by ma- 
jor flood basalt eruptions and massive faunal depletions 
[la ]. The authors suggested that large thermal plumes 
ascending through the mantle favor reversals and sub- 
sequently produce large eruptions. In the light of our 
work, it is tempting to associate to the thermal plumes 
(which provide a localized thermal forcing at the sur- 
face of the CMB) with an enhanced deviation from the 
flow equatorial symmetry, which results according to the 
above description in an increase of reversal frequency and 
therefore ends superchrons. 

In contrast, another scenario has been proposed in 
which the Earth's dipole is coupled to an octupole, i.e., 
another mode with a dipolar symmetry [ll|. This does 
not require additional constraint on the flow in the core 
in the framework of our model. In any case, the existence 
of two coupled modes allows the system to evolve along 
a path that avoids B = 0. In physical space, this means 
that the total magnetic field does not vanish during a 
reversal but that its spatial structure changes. 

Numerical simulations of MHD equations [l|| or of 
mean field models have displayed reversals that seem to 
involve "transitions between the steady and the oscilla- 
tory branch of the same eigenmode" [13. That situa- 
tion can be obtained in the vicinity of a codimension-two 
bifurcation with a double zero eigenvalue and only one 
eigenmode. This type of bifurcation also exists in our 
model ([2]) but requires tuning of two parameters /i and 
v. This is not necessary for the scenario of reversals we 
have described. Other features of reversals observed in 
numerical simulations at magnetic Prandtl number of or- 
der one, such as mechanisms of advection/amplification 
of the field due to localized flow processes are not 
described by our model which requires the limit of small 
magnetic Prandtl number (relevant to the Earth's core). 

Equations similar to ([3]) have been studied in a variety 



of problems, for instance for the orientation of a rigid 
rotator subject to a torque [Hj], used as a toy model for 
the toroidal and the poloidal field of a single dynamo 
mode. Indeed, symmetries constrain the form of the 
equation for 9 even though the modes and the physics 
involved are different. We emphasize that the above 
scenario is generic and not restricted to the equation 
considered here. Limit cycles generated by saddle-node 
bifurcations that result from the coupling between two 
modes occur in Rayleigh-Benard convection [2(J. A 
similar mechanism can explain reversals of the large scale 
flow generated over a turbulent background in thermal 
convection or in periodically driven flows [2l[. We 
have proposed a scenario for reversals of the magnetic 
field generated by dynamo action that is based on the 
same type of bifurcation structure in the presence of 
noise. It offers a simple and unified explanation for 
many intriguing features of the Earth's magnetic field. 
The most significant output is that the mechanism 
predicts specific characteristics of the field obtained 
from paleomagnetic records of reversals and from recent 
experimental results. Other characteristic features such 
as excursions as well as the existence of superchrons 
are understood in the same framework. Below the 
threshold of the saddle-node bifurcation, fluctuations 
drive random reversals by excitability. We also point out 
that above its threshold, the solution is roughly periodic. 
It is tempting to link this regime to the evolution of 
the large scale dipolar field of the Sun (which reverses 
polarity roughly every 11 years). Recent measurements 
of the Sun surface magnetic field have shown that two 
components oscillate in phase-quadrature [22| . This 
would be coherent with the oscillatory regime above the 
onset of the saddle-node bifurcation if these components 
correspond to two different modes. 

We thank our colleagues from the VKS team with 
whom the data published in [3[ have been obtained. 
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